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Abstract 

We investigate the properties of astrophysical electromagnetic cascades in matter, 
photon gas and magnetic fields, and discuss similarities and differences between 
characteristics of electron-photon showers developed in these 3 substances. We ap- 
ply the same computational technique based on solution of the adjoint cascade 
equations to all 3 types of cascades, and present precise numerical calculations of 
cascade curves and broad-band energy spectra of secondary electrons and photons 
at different penetration lengths. 
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1 Introduction 

Relativistic electrons - directly accelerated, or being secondary products of 
various hadronic processes - may result in copious 7-ray production caused 
by interactions with ambient targets in forms of gas (plasma), radiation and 
magnetic fields. In different astrophysical environments 7-ray production may 
proceed with high efficiency through bremsstrahlung, inverse Compton scat- 
tering and synchrotron (and/or curvature) radiation, respectively. Generally, 
7-ray production is effective when the cooling time that characterizes the rate 
of the process does not significantly exceed (i) the source (accelerator) age (ii) 
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the characteristic time of non-radiative losses caused by adiabatic expansion 
of the source or particle escape and (Hi) the cooling time of competing radia- 
tion mechanisms that result in low-energy photons outside the 7-ray domain. 
As long as the charged particles are effectively confined in the 7-ray produc- 
tion region, at some circumstances these condition could be fulfilled even in 
environments with a relatively low gas and photon densities or weak mag- 
netic field. More specifically, the 7-ray production efficiency could be close 
to 1 even when t rad R/c (R is the characteristic linear size of the pro- 
duction region, c is the speed of light). In such cases the secondary 7-rays 
escape the source without significant internal absorption. Each of the above 
mentioned 7-ray production mechanisms has its major "counterpart" - 7-ray 
absorption mechanism of same electromagnetic origin resulting in electron- 
positron pair production in matter (the counterpart of bremsstrahlung) , in 
photon gas (the counterpart of inverse Compton scattering), and in magnetic 
field (the counterpart of synchrotron radiation). The 7-ray production mecha- 
nisms and their absorption counterparts have similar cross-sections, therefore 
the condition for radiation t ra d > R/c generally implies small optical depth for 
the corresponding 7-ray absorption mechanism, r abs < 1. But in many astro- 
physical scenarios, in particular in compact galactic and extragalactic objects 
with favorable conditions for particle acceleration, the radiation processes pro- 
ceed so fast that t ra d > R/c. At these conditions the internal 7-ray absorption 
becomes unavoidable. If the 7-ray production and absorption processes oc- 
cur in relativistic regime, namely when (i) E lfi > 10 3 m e c 2 in the hydrogen 
gas, (ii) E lfi e ^> m 2 c 4 in photon gas (often called Klein-Nishina regime; e is 
the average energy of the target photons), or (iii) (E lfi / m c c 2 )(H / H CT ) 3> 1 
in the magnetic field (often called quantum regime; H cr ~ 4.4 x 10 13 G is 
the so-called critical strength of the magnetic field), the problem cannot be 
reduced to a simple absorption effect. In this regime, the secondary electrons 
produce new generation of high energy 7-rays, these photons again produce 
electron-positron pairs, so the electromagnetic cascade develops. 

The characteristics of electromagnetic cascades in matter have been compre- 
hensively studied, basically in the context of interactions of cosmic rays with 
the Earth's atmosphere (see e.g. [1]) as well as for calculations of performance 
of detectors of high energy particles (e.g. [2]). The theory of electromagnetic 
cascades in matter can be applied to some sources of high energy cosmic 
radiation, in particular to the "hidden source" scenarios like massive black 
holes in centers of AGN or young pulsars inside the dense shells of recent 
supernovae explosions (see e.g. [3]). Within another, so-called "beam dump" 
models (see e.g. Ref. [3,5]) applied to X-ray binaries, protons accelerated by 
the compact object (a neutron star or black hole), hit the atmosphere of the 
normal companion star, and thus result in production of high energy neu- 
trinos and 7-rays [4]. In such objects, the thickness of the surrounding gas 
can significantly exceed 100 g/cm 2 , thus the protons produced in the central 
source would initiate (through production of high energy 7-rays and electrons) 
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electromagnetic showers. These sources perhaps represent the "best hope" of 
neutrino astronomy, but they are generally considered as less attractive targets 
for gamma-ray astronomy. However, the 7-ray emission in these objects is not 
fully suppressed. The recycled radiation with spectral features determined by 
the thickness ("grammage") of the gas shell, should be seen in 7-rays in any 
case, unless the synchrotron radiation of secondary electrons dominates over 
the bremsstrahlung losses and channels the main fraction of the nonthermal 
energy into the sub gamma-ray domains. 

The development of electromagnetic cascades in photon gas and magnetic 
fields is a more common phenomenon in astrophysics. In photon fields such 
cascades can be created on almost all astronomical scales, from compact ob- 
jects like accreting black holes [6-10], fireballs in gamma-ray bursts [11] and 
sub-pc jets of blazars [12] to large-scale (up to > 100 kpc) AGN jets [13,14] and 
> 1 Mpc size clusters of galaxies [15]. Electromagnetic cascades in the inter- 
galactic medium lead to formation of huge (> 10 Mpc) nonthermal structures 
like hypothetical electron-positron pair halos [16]. Finally, there is little doubt 
that the entire Universe is a scene of continuous creation and development of 
electromagnetic cascades. All 7-rays of energy > a few GeV emitted by astro- 
physical objects have a similar fate. Sooner or later they terminate on Hubble 
scales due to interactions with the diffuse extragalactic background. Since the 
energy density of 2.7 K CMBR significantly exceeds the energy density of 
intergalactic magnetic fields, these interactions initiate electron-photon cas- 
cades [17]. The superposition of contributions of 7-rays from these cascades 
should constitute a significant fraction of the observed diffuse extragalactic 
background. 

Bonometto and Rees [18] perhaps where the first who realized the astrophys- 
ical importance of development of electron-photon cascades supported by 7-7 
pair-production and inverse Compton scattering in dense photon fields. When 
the so called compactness parameter [19] / = (L/R)(<jT</m e c 3 ) (L is the lu- 
minosity and R is the radius of the source) is less than 10, then the cascade 
is developed in the linear regime, i.e. when the soft radiation produced by 
cascade electrons do not have a significant feedback effect on the cascade 
development. In many cases, including the cascade development in compact 
objects, this approximation works quite well. 

The first quantitative study of characteristic of linear cascades in photon fields 
has been performed using the method of Monte Carlo simulations [20]. Gen- 
erally, the kinetic equations for the cascade particles can be solved only nu- 
merically. However, with some simplifications it is possible to derive useful 
analytical approximations [7,8] which help to understand the features of the 
steady-state solutions for cascades in photon fields. 

The cascade development in the magnetic field is a key element to understand 
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the physics of pulsar magnetospheres [21,22], therefore it is generally treated as 
a process associated with very strong magnetic fields. However, such cascades 
could be triggered in many other (at first glance unusual) sites like the Earth's 
geomagnetic field [23-25], accretion disks of massive black holes [26], etc. In 
general, the pair cascades in magnetic fields are effective when the product of 
the particle (photon or electron) energy and the strength of the field becomes 
close to the "quantum threshold" of about H cr i t m e c 2 ~2x 10 7 TeV- G, unless 
we assume a specific, regular field configuration. An approximate approach, 
similar to the so-called approximation A of cascade development in matter [1], 
has been recently applied by Akhiezer et al. [27]. Although this theory quite 
satisfactorily describes the basic features of photon-electron showers, it does 
not provide an adequate accuracy for quantitative description of the cascade 
characteristics [24]. Note that in both studies approximate cross-sections for 
magnetic bremsstrahlung and magnetic pair-production have been used, thus 
reducing the validity of the results to the limit EH » 10 7 TeV • G (E is the 
minimum energy of secondary particles being under consideration). 

As long as we are interested in the one-dimensional cascade development 
(which seems to be quite sufficient for many astrophysical purposes), all 3 
types of cascades can be described by the same integro-differential equations 
like the ones derived by Landau and Rumer [28], but in each case specifying 
the cross-sections of relevant interaction processes. The solution of these equa- 
tions in a broad range of energies is however not a trivial task. In this paper 
we present the results of our recent study of cascade characteristics in 3 sub- 
stances - matter, photon gas and magnetic field - with emphasis on the analysis 
of similarities and differences between these 3 types of cascades. For quanti- 
tative studies of these characteristics we have chosen the so-called technique 
of adjoint cascade equations. Although this work has been initially motivated 
by methodological and pedagogical objectives, some results are rather original 
and may present practical interest in certain areas of high energy astrophysics. 



2 Technique of adjoint cascade equations 

The results of this study are based on numerical solutions of the so-called 
adjoint cascade equations. The potential of this computational technique, in 
particular its possible applications to different problems of cosmic ray physics 
has been comprehensively described in Ref.[29]. Below we briefly discuss the 
main features of solution of the adjoint cascade equations. For our purposes it 
is quite sufficient to consider only the longitudinal development of electromag- 
netic cascades neglecting the emission angles of secondary particles and, thus, 
assuming that all cascade particles are moved along the shower axis. Besides, 
we neglect the differences in interactions of electrons and positrons, i.e. do not 
consider positron annihilation, both in flight and after thermalization in the 
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ambient plasma. This is a quite good approximation as long as we are inter- 
ested in 7-ray energies exceeding a^ 1 m e c 2 ~ 70 MeV (aif = 1/137 is the fine 
structure constant). Under these assumptions, the system of adjoint cascade 
equations reads 

% + f/ X e~ J W ee (e,e')f(t,e')de'- j W ej (e, e')g(t, e')de' = F(t, e),(l) 
JL + g/X^- J W ie (e,e')f(t,e')de'- J W J7 (e, e')g(t, e')de' = G(t, e).(2) 



The adjoint functions fit, e) and g(t, e) in Eqs. (1) and (2) describe the contri- 
butions (averaged over random realizations of the cascade development) from 
cascade initiated by a primary electron (/) or a photon (g) of energy e (all 
energies are expressed in units m e c 2 ); t is the cascade penetration depth; e t h 
is the minimum energy of cascade particles being under consideration. 

The parameter \ a in Eqs. (1) and (2) defines the mean free path length of 
cascade particles of type a (a — e or 7). It can be expressed through the 
total cross-sections (cr a ) of interaction with the ambient medium. In the most 
general case, when the medium consists of matter (M), magnetic field (F) and 
photon gas (G), 



X a — 



n { M) a^+a^+nPaP]- L (3) 



where and are the number density of the matter atoms and the 

ambient photons, respectively. The parameters 

Wctie, e') = ni M) wi M p \e, e>) + w%\e, e>) + n^w^e, e') (4) 



are differential cross-sections over the energy e' of the secondary particle of 
type P {(5 — e or 7). They are normalized as 

/ W a ,(e, e')de' = r^n™*™ + n^a^ + n<$*S» (5) 



where n a/3 , n a/3 and n a/3 are the mean multiplicities of secondary particles 
of type P produced by a particle of type a. 

Properties of the particle "detector" located at the depth t are defined by the 
right hand side functions F and G in Eqs. (1) and (2) and by the boundary 
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conditions f(t = 0,e) and g(t = 0,e). For example, if the "detector" measures 
the number of cascade electrons above some threshold energy, e > e t h, then 

F(t,e) = G(t,e) = 0; f(t = 0, e) = H(e - e th ), g(t = 0,e) = (6) 

where H(x) = 1 for x > and H(x) = for x < 0. Analogously, if the 
"detector" measures the number of cascade photons with energy e > e t h, then 

F(t, e) = G(t, e) = 0; f(t = 0, e) = 0, g(t = 0, e) = H(e - e th ). (7) 



For solution of the adjoint cascade equations we use the numerical method 
proposed in Ref. [30]. To apply this approach to Eqs. (1) and (2) we introduce 
an increasing subsequence {e^} = e = e t h, £i, ...,£&,... of energy points, and 
corresponding values of adjoint functions {fk(t),9k(t)}- Let's write now the 
Lagrange polynomial interpolation formulae for an approximate presentation 
of functions / and g inside the energy intervals Ae^ = (£fc_i, £&) for k — 1, 2, . . . 
through their values at the support points {£&}, 

f(e,t)^f(e,t)= £ fiW^e), g(e,t) ~ g(e,t) = E ^(t)^.(e); 

j=k~n j=k—n 
n 

L kj( £ ) = T[( £ - £ k-r)(£j -£fe-r) _1 , r^k-j, eeAe k (8) 

r=0 



where n = n(k) = min(A^, k) and iV is the maximum available power of the 
Lagrange polynomial. 

For the first k energy intervals Ae fe the support values from fo,go to fk,gk 
are used; in this case the polynomial power n is equal to k. For intervals with 
k > N the support values f k „ N ,g k - N , fk-N+i, 9k-N+i, ■ ■ ■, fk,9k are involved 
in the interpolation with use of the power n = N. 

Let's adopt now e = e k in Eqs.(l) and (2), and present the integral members 
of these equations in the form of the sum over the energy intervals Ae, (i = 
1,2,..., k): 

/ W a p(e k ,E)---de = J2 W al3 (e k ,e)---de, a,f3 = e,^. (9) 

i— 1 



Then, after some simple calculations, we find the following equations 

^ + A k f k - B k g k — F' k , ^ - C k f k + D k g k = G' k (10) 
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where 



A k — 1/X e (e k ) — a kk , B k — b kk , Ck — Ckk, D k — 1/A 7 (£jfc) — d kk ,(ll) 

k-l k-l 

F 'k = F k + Y, \ a 3kfi + b jk gj] , G' k = G k + Y [cjkfj + d jk gj] . (12) 

j=0 i=o 



The coefficients a^, frj/t, Cjfe and rfj^ can be expressed through cross-sections of 
relevant processes: 

i=lg j _ 1 s=i—n 
k e j i 

b jk = Y J W ei (e k ,e)Ll(e) £ , etc. (13) 



i=l F . 



with <5j S = 1 if j = s, and Sj S = if j ^ s. 

Thus, this method allows us to reduce the integro-differential Eqs. (1) and (2) 
to the system of ordinary differential equations (10). The solution of Eq. (10) 
can be obtained in two steps: 

(1) for e = eo = e t h in formulae (1) and (2), we find the following equations 
for functions f (t) and g (t) 

^ + fo(t)/\ e (e ) = F(t,s ), ^ + 9o (t)/\ 1 (e ) = G(t,s ) (14) 



(2) after solving Eqs. (14) we can calculate the functions /i (t) , <7i (t) , because 
for k — 1 Eq. (10) contains only fo(t),go(t), fi(t) and <?i(t); after that one can 
find by the same way the functions f2(t),g2(t), f3(t),g3(t), etc. 

For a fixed value of k, we introduce an increasing subsequence of the depth 
values {*,} (t = 0,tj = tj-i+7}) and correspondingly f k (t t ) = f kh g k (U) = g Kl . 
For each interval (t;_i,tj) one can evaluate fk,i,g k ,i by solving Eq. (10) for 
which the values f k ,i-i, g k ,i-i serve as boundary conditions. This leads to 

i 

fk,i = (A fc - Aife) -1 J2(- l Y{ ex P( X vkTi)[(D k + \ v k)fk,i-i + B k g k ^i]+ 

ti 

+ J ew[Kk(ti-T)][(D k + \ vk )F^T)+B k G' k (T)]}dT, (15) 
ti-i 



7 



1 

9k,i = (Aofc - Aife) -1 ^2(-l) u {exp(\ vk Ti)[(A k + \ vk )gk,i-\ + C k fk,i-i] + 

ti 

+ J eMKk{ti-r)}[{A k + \ uk )G' k {r) + C k Fl{r)}}dT (16) 

ti-i 

where 

Kk = + D k ) + (-iy[(A k - D k f + AB k C k ]^ 2 }. (17) 

Note that for the homogeneous environment these relations are valid for an 
arbitrary value of 77, otherwise they can be applied only for 77 <C U. 

The knowledge of f k ,o and g k ,o (these quantities can be calculated on the basis 
of boundary conditions like Eq. (6) or Eq. (7) ) and the multiple application 
of Eqs (15) - (17) allow to calculate quantities f k ,i,gk,i f° r I — 1,2, etc. 

The approach of solution of adjoint cascade equations described above provides 
an accuracy better than a few per cent and gives results for an arbitrary region 
(£th, ^max) of the primary energy in just one run of calculations. Also we notice 
that the computational time consumed by this method does not exceed on 
average a few minutes for a 1 GHz PC type computer and increases only 
weakly (logarithmically) with the primary energy. In summary, the important 
feature of this technique is its flexibility to describe the cascade processes 
in 3 different substances. It allows large number of calculations with a good 
accuracy throughout very large energy intervals of both primary and secondary 
energies. This is an important condition for the quantitative analysis and for 
clear understanding of similarities and differences in cascade development in 
environments dominated by matter, photon gas or magnetic field. 



3 Elementary processes 

The above described technique requires specification of elementary processes 
that initiate and support the cascade development. If we are interested in 
the longitudinal development of cascades, the input should consist of total 
cross-sections of interactions and the differential cross-sections as functions 
of energy but integrated over the emission angles of secondary particles. In 
many astrophysical situations, especially at very high energies, this is a fair 
approximation, given the very small (of order ~ m c c 2 /E 1) emission angles 
of secondary products. The one-dimensional treatment of the cascade develop- 
ment perfectly works, if electrons move along the lines of the regular magnetic 



8 



field. Otherwise, the deviation of particles from the cascade core is determined 
by reflection of secondary electrons by magnetic inhomogeneities rather than 
the emission angles. In these cases the diffusion effects must be appropriately 
incorporated into the cascade equations. This question is beyond the present 
paper. 

All processes involved in the electromagnetic cascades are well known. The 
cross-sections of these processes have been calculated and comprehensively 
studied with very high precision within the quantum electrodynamics. 

In the case of cascades in matter we take into account the ionization losses 
and bremsstrahlung for electrons (positrons), and the pair production, Comp- 
ton scattering and photoelectric absorption for photons. At extremely high 
energies the so-called Landau- Pomeranchuk-Migdal (LPM) effect that results 
in suppression of the bremsstrahlung and pair-production cross-sections, may 
have significant impact on the cascade development. The energy region where 
the LPM effect becomes noticeable, depends on the atomic number of the am- 
bient matter. In the hydrogen dominated medium this energy region appears 
well beyond E > 10 20 eV, therefore in many astrophysical scenarios the LPM 
effect can be safely ignored. 

Two pairs of coupled processes - the inverse Compton scattering and 7-7 pair 
production in the photon gas and the magnetic bremsstrahlung (synchrotron 
radiation) and magnetic pair production in the magnetic field - determine 
the basic features of cascade produced in radiation and magnetic fields. At 
extremely high energies the higher order QED processes may compete with 
these basic channels. Namely, when the product of energies of colliding cascade 
particles (electrons or photons) E and the background photons e significantly 
exceed 10 5 — 10 6 m 2 c 4 , the processes 77 — * e + e~e + e~[31] and — > eje + e~ 
[32] dominate over the single (e + , e~) pair production and the Compton scat- 
tering, respectively. For example, in the 2.7 K CMBR the first process stops 
the linear increase of the mean free path of highest energy 7-rays around 
10 21 eV, and puts a robust limit on the mean free path of 7-rays of about 100 
Mpc. Analogously, above 10 20 eV the second process becomes more important 
than the conventional inverse Compton scattering. Because 77 — > 2e + 2e~ and 
e7 — > e / -/e + e~ channels result in production of 2 additional electrons, they 
significantly change the character of the pure Klein-Nishina pair cascades. 
However, an effective realization of these processes is possible only at very 
specific conditions with an extremely low magnetic field and narrow energy 
distribution of the background photons. 

Because this study has methodological objectives, below we do not include 
these processes in calculations of cascade characteristics. This allows us to 
avoid unnecessary complications and make the analysis more transparent. For 
the same reason we do not include in this study the effect of the photon 
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splitting [33] which becomes quite important in pulsars with magnetic field 
close to # crit = 4.4 x 10 13 G (see e.g. Ref. [22]). 



3.1 Total cross-sections 



In Figs. 1,2 and 3 we show the photon- and pair-production cross-sections in 
hydrogen gas and in radiation and magnetic fields, respectively. The energies 
of electrons and 7-rays are expressed in units of m e c 2 . 




Fig. 1. Total cross-sections of the bremsstrahlung (o"b r ) and pair production (<r p ) pro- 
cesses in hydrogen gas normalized to the asymptotic value (er ) of the pair produc- 
tion cross-section at £0 —> 00 • The bremsstrahlung cross-sections are calculated for 
secondary 7-rays produced with energies exceeding (1) the pair-production thresh- 
old, £th = 2; (2) the critical energy, e t h — 700; (3) half of the energy of the primary 
electron, e t h = eo/2. 

In Fig. 1 the total cross-sections for matter are presented. They are normalized 
to the asymptotic value of the pair production cross-section at Eo — > 00: 

a = 7/9 x 4 ai r 2 c Z(Z + 1) ln ( 183Z " 1/3 ) (18) 
° 1 c v ; 1 + 0.12 Z/82 2 v ; 



where Z is the charge of the target nucleus, r e is the classical electron radius. 
This actually implies introduction of the so-called radiation length 



X ( M) = 7/9 



n. 



(M) 1 
°o 



(19) 
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Table 1 

The asymptotic behaviour of the total cross-sections of electron and photon inter- 
actions at high energies. 



Environment 


Matter 


Photon gas 


Magnetic field 




~ lneo 


~ (no ■ lnKo)^ 1 


-1/3 


(7-y 


~ const 


~ (k • hi/to) -1 


-1/3 

~Xo 



Xq has a meaning of the average distance over which the ultrarelativis- 
tic electron loses all but 1/e of its energy due to bremsstrahlung. The same 
parameter approximately corresponds to the mean free path of 7-rays. There- 
fore, the cascade effectively develops at depths exceeding the radiation length. 
Usually the radiation length is expressed in units g/cm 2 . For the hydrogen 
gas Xq M ^ = 63 g/cm 2 . The second important parameter that characterizes 
the cascade development is the so-called critical energy below which the ion- 
ization energy losses dominate over bremsstrahlung losses. In the hydrogen gas 
e CT ~ 700. Effective multiplication of particles due to the cascade processes is 
possible only at energies e > e CT . At lower energies electrons dissipate their 
energy by ionization rather than producing more high energy 7-rays which 
would support further development of the electron-photon shower. 

The bremsstrahlung differential cross-section has a l/e 7 type singularity at 
e 7 — > (e 7 is the energy of the emitted photon), but because of the hard 
spectrum of bremsstrahlung photons the energy losses of electrons are con- 
tributed mainly by emission of high energy 7-rays. In Fig. 1 we show the 
bremsstrahlung total cross-sections calculated for 3 different values of mini- 
mum energy of emitted photons: e t h = 2, e cr and e c /2. The first value cor- 
responds to the cross-section of production of all 7-rays capable to produce 
electron-positron pairs. The second value corresponds to the cross-section of 
production of 7-rays produced above the critical energy, and thus capable to 
support the cascade. And finally, the third value corresponds to the cross- 
section of production of the "most important" 7-rays which play the major 
role in the cascade development. It is seen from Fig. 1 that while for £ t h = 2 
the total cross-section of pair-production is an order of magnitude lower com- 
pared to the bremsstrahlung cross-section, for £ th = e c /2 the cross-sections of 
two processes become almost identical at energies e > 100. 

Both cross-sections grow significantly with energy until £0 ~ 10 3 . At higher 
energies the pair-production cross-section becomes energy-independent, but 
for £ th = 2 the bremsstrahlung cross-section continues to grow, although very 
slowly (logarithmically). 

In Fig. 2 we show the total cross-sections of the inverse Compton scatter- 
ing and pair-production for the mono-energetic and Planckian isotropic pho- 
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Fig. 2. Total cross-sections of the inverse Compton scattering and pair production. 
Two spectral distributions for the ambient photon gas are assumed - mono-energetic 
with energy c^o (curves 1 and 3) and Planckian with the same mean photon energy 
ujq ~ 3kT/m c c 2 . 

ton fields, normalized to the Thompson cross section <7t = 8/37rr 2 = 6.65 x 
1CT 25 cm 2 . Both cross sections depend on the product k = e u of the pri- 
mary (eq) and ambient (u = e/m c c 2 ) photon energies. At k — > the inverse 
Compton cross-section a\c — > cr T . At high energies it decreases with k , as 
oc (/to • In/to) -1 . 

In the mono-energetic radiation field the pair production has a strict kinematic 
threshold at k = 1. At small values of k < 2 the cross-section rapidly 
increases with k achieving its maximum of about ~ 0.2<r T at k ~ 3.5 — 4, and 
decreases with further increase of Kq as oc (/t Tn/t ) _1 . Thus in this regime the 
pair-production cross-section behaves quite similar to the Compton scattering 
in the Klein-Nishina regime, but its absolute value is higher by a factor of 1.5. 

In the case of interactions of electrons and photons with the magnetic field 
it is more convenient to introduce, instead of standard total cross-sections, 
the interaction probabilities [24]. But in the literature this parameter is still 
formally called as cross-section. These probabilities normalized to the strength 
of the magnetic field are shown in Fig. 3. The probabilities of both processes 
- photon production by electrons (synchrotron radiation) and pair production 
by photons - depend only on the parameter Xo — £oH/H cr , where H is the 
component of magnetic field intensity perpendicular to the velocity of the 
particle. The parameter xo = H/ H ct Sq is an apparent analog of the parameter 
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10"' 1 10 io 2 io 3 io 4 io 5 

Zo=H/H cr e 

Fig. 3. Total cross-sections of the synchrotron radiation and magnetic pair produc- 
tion. Solid curves - our calculations; points - calculations of Anguelov et.al. [24]. 



Table 2 

The singularities in the differential cross-sections of cascade processes. 



Elementary process Singularity point Behaviour 



Bremsstrahlung 




= 




< 


Inverse Compton scattering 




= £o 


~ (e 




Pair production in photon gas 


e e 


= 






Pair production in photon gas 


e e 


= £o 




-Se)- 2 


Synchrotron radiation 


By 


= 


~ i 


-2/3 

'7 



kq. While the probability of the synchrotron radiation at Xo ^ 1 is constant, 
the probability of the pair production drops dramatically below xo — 1- After 
achieving its maximum at Xo — 10, the probability of the pair production 
decreases with x as oc Xq 1 ^ ■ The same behaviour at large x nas a l so the 
probability of synchrotron radiation, but the absolute value of the latter always 
exceeds by a factor of 3 the probability of the pair production. 
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Table 3 

The mean portion of primary energy (e 7 /eo) transferred to the secondary photon in 
the inverse Compton scattering (ics) and the synchrotron radiation (syn) processes 
via the values of parameters kq = Eqluo and %o = £qH/H ct , respectively. 



«0, XO 


0.01 


0.1 


1 


10 2 


10 4 


10 6 


(e 7 /eo)ics 


0.014 


0.099 


0.358 


0.760 


0.867 


0.910 


(e 7 /£ )syn 


0.44 • 10~ 2 


0.033 


0.118 


0.241 


0.250 


0.250 



Table 4 

The mean portion of primary energy (e e /eo) transferred to the leading secondary 
particle for the processes of pair production in a mono-energetic photon gas (G) and 
in the magnetic field (F) via the values of parameters kq = eo^o and xo = £oH/H CT , 
respectively. 



«o, Xo 


1 


3 


10 


10 2 


10 4 


10 6 


(e e /e )G 


0.500 


0.701 


0.797 


0.891 


0.948 


0.966 


(e e /e )F 


0.634 


0.693 


0.746 


0.782 


0.824 


0.825 



3.2 Differential cross-sections of cascade processes 

The differential cross-sections of cascade processes are presented in Figs. 4,5 
and 6. The bremsstrahlung and pair-production cross-sections are from Ref. 
[34,35]. The cross-sections for the inverse Compton scattering and pair pro- 
duction in the mono-energetic isotropic photon field are from Ref. [35] and 
Ref. [36], respectively. The cross- sect ions of processes in the magnetic field are 
from Ref. [27]. 

All three cross-sections of pair-production are symmetric functions in respect 
to the point x = e e /eo = 0.5. The photon production processes are asymmetric 
functions. The bremsstrahlung, synchrotron radiation, the inverse Compton 
scattering and the pair production in the photon gas have singularities. The 
location of singularities and the cross-section behaviour near the singularity 
points are summarized in Table 2. 

The character of cascade development is largely determined by the fraction 
of energy of electrons and photons lost per interaction. In Table 3 we present 
mean fractions of the electron energy transferred to the secondary photons in 
the inverse Compton and synchrotron radiation processes at different values of 
and Xo- 111 the classical regime (kq <C 1 and xo <S 1) the mean energy lost 
by the electron per interaction is very small, Ae/e <C 1. In the "quantum" 
regime (kq ^> 1 and Xo ^ 1) the interactions have a catastrophic charac- 
ter; the secondary photons get a significant fraction of the energy of parent 
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Fig. 4. Differential cross-sections of the bremsstrahlung (upper panel) and pair pro- 
duction (bottom panel) processes in hydrogen. The cross-sections are normalized to 
one radiation length. The energies of primary electrons and 7-rays Eq are indicated 
at the curves. 

electrons. In the photon gas at k,q ^> 1 this fraction exceed 0.5, approach- 
ing asymptotically to 1. In the magnetic field the energy transfer is smaller; 
at xo ~ 1 it is approximately 0.1, and asymptotically approaches to 1/4 at 
extremely large xo- 

The pair production processes in all 3 substances have (by definition) catas- 
trophic character (the photon always disappears). Since the differential spectra 
of secondary electrons are quite flat with increase towards the maximum en- 
ergy e e — > e 1 (in the ultrarelativistic regime), these processes proceed with 
formation of leading electron with energy x = e e /e 7 — * 1 in the photon gas 
and somewhat smaller (~ 0.8) in the magnetic field (see Table 4). 
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Fig. 5. Differential cross-sections of the inverse Compton (upper panel) and 
pair production (bottom panel) processes for the case of a mono-energetic gas 
of ambient photons. The parameters e imax , £ e min and e cmax are defined as 
e 7 max = 4eo(^o/l + 4k ) and £ e min,emax = 0.5eo(l T v 7 ! ~ V K o)- Different values of 
the parameter kq = £qujq are indicated at the curves. 



4 Cascades 



In this section we discuss and compare the so-called cascade curves and the en- 
ergy spectra of electrons and photons for showers produced in matter, photon 
gas and magnetic field. 
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Fig. 6. Differential cross-sections of synchrotron radiation (upper panel) and the 
magnetic pair production (bottom panel) normalized to the total cross-sections of 
these processes. Different values of parameter xo = H/H cr eo are indicated at the 
curves. 

4-1 Cascade curves 



The cascade curve Np(t, e , £ t h) describes the average number of cascade elec- 
trons ((3 = e) or photons {(5 = 7) above e t h, as a function of the penetration 
depth t. 



4.1.1 Matter 

In Fig. 7 we show the cascade curves for electrons calculated for different values 
of the primary (e ) and threshold (e t h) energies using the adjoint equation 
technique. For comparison we also show the cascade curves calculated within 
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the approximation A of the analytical theory of electromagnetic cascades [1] 
(valid for e t h ^ £ CT ), as well as the cascade curves obtained by Monte Carlo 
simulations using the ALTAI code[37]. The results obtained by 3 different 
methods are in a good agreement with each other. 

One of the basic parameters characterizing the cascade development is the 
depth t m (expressed in units of radiation length Xq M ^) at which the cascade 
curve achieves its maximum. Generally t m grows logarithmically with the pri- 
mary energy. In particular, for light materials (with the nucleus charge number 
Z < 10) and for e £cr, the parameter t m for electrons of a cascade initiated 
by a primary photon can be approximated as [1]: 

t c m ~ ln(e /eth) for £th > £ cr ; t m ~ ln(e /e cr ) for e th -> 0. (20) 

The primary electron interacts with matter somewhat earlier than the primary 
7-ray (see Fig. 1), therefore + 0.7. 

Another important parameter is the total number of electrons at the shower 
maximum: N max = N e (t m , Eq, £ t h)- This number is approximately proportional 
to the primary energy. In the energy region £ t h ^> s CT the electron number 
Acm ~ ^th. 1 ; however it does not depend on e t h, if £th — *■ 0. 

In Fig. 8 we present the ratio of the number of the cascade photons to electrons 
as a function of the penetration depth. For e th e cr , the number of cascade 
photons A 7 is comparable to the electron number N e . At small threshold 
energies, (e t h ^ ^cr) the number of 7-rays considerably exceeds the number of 
electrons. This is explained by the break of symmetry between the electrons 
and 7-rays caused by ionization losses of electrons. 

4-1.2 Photon gas 

For description of the cascade development in the mono-energetic photon gas it 
is convenient to introduce, analogously to the cascade in matter, the radiation 
length Aq G ^ in the following form 

1 ^0 , (21) 

where n 1 ^' is the number density of photons. Apparently, corresponds, 
with accuracy up to logarithmic terms, to the mean free path of 7-rays in the 
photon field at k 3> 1. 

Note that unlike to the radiation length in matter, the primary energy en- 

(G) 

ters, through the parameter k = Equjo, in Xq . Obviously, for the ambient 



y G) _ 
A - 
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Fig. 7. Cascade curves of electrons for showers initiated by primary electrons (solid 
curves) and photons (dashed curves) . The calculations are performed for the follow- 
ing primary energies Eq = 2 • 10 8 (curve 1), 2 • 10 7 (curve 2), 2 • 10 4 (curve 3) and 
the ratio e t h/£ C r = 125 (curves 1 and 2), 0.05 (curve 3). For comparison, the results 
derived from the analytical cascade theory [1] (boxes) and by simulations with the 
ALTAI code [37] (triangles) are also shown 

photon gas with a broad band energy distribution this parameter becomes 
meaningless. At the same time in a narrow band radiation fields, e.g. with 
Planckian distribution, this parameter can work effectively by substituting 
Cu w 3kT/m c c 2 . Note that for the black-body radiation the density of pho- 
tons is determined by temperature, therefore 

X^ BB) = 3/7ralr (kT/m c c 2 y 2 E ~ 6.9 x lCT 7 \kT/m c c 2 y 2 e cm (22) 



In particular, in the 2.7 K CMBR X^ a> ~ 3.3 x 10 12 e cm. 

The results shown in Fig. 9 are cascade curves of electrons and photons in 
the blackbody photon gas calculated for the fixed value of K th = e^uJo = 1 
corresponding to the threshold energy of cascade particles, £ t h = m c c 2 /3kT, 
but for different values of the parameter kq, i.e. for different primary photons 
energies e = k ■ m c c 2 /3kT. 

We can see that the depth t m of the shower maximum measured in units of 
radiation length Xq depends rather weakly on the primary energy £ - It 
ranges within ~ 1 -f- 2. This means that in geometrical units of length t m is 
approximately a linear function of Eq. This is explained by the approximately 
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Fig. 8. The depth dependence of the N^/N e ratio for cascades initiated by primary 
electrons (solid curves) and photons (dashed curves). Primary energy eq = 2 ■ 10 8 ; 
£th/£cr = (curves 1), 125 (curves 2). 

linear decrease of the cross- sect ions and correspondingly by the increase of 
mean free path of electrons and 7-rays interacting with the ambient photons 
in the Klein-Nishina (quantum) regime (see Fig. 2). 

The number of cascade particles in the Klein-Nishina regime increases with 
Eq slowly. Even near the shower maximum this number does not exceed a few 
particles. This is explained by an extremely high efficiency of conversion of en- 
ergy of the photon to the leading electron at 7-7 interactions, and vice versa - 
the electron energy to the upscattered photon at the Compton scattering (see 
Figure 5 and Tables 3,4). As a result, the energy of the second (secondary) 
particle appears too small for noticeable contribution to the cascade develop- 
ment. Since in this regime the cross-sections of the Compton scattering and 
photon-photon pair production are quite similar, the cascade process in the 
zeroth approximation can be considered as propagation of a single composite 
7/e" particle which spends 2/5 of its time in the "7 state" and 3/5 in the "elec- 
tron state" (these times are determined by the ratio 1.5 of the corresponding 
cross-sections in the Klein-Nishina limit). 

Fig. 10 illustrates dependence of the cascade curve on the threshold energy e th 
with the corresponding parameter fc t h below and above 1. With reduction of 
«4h, the cascade curves of both 7-rays and electrons increase, especially in the 
regime of fc t h < 1- This has a simple explanation. Although in this (Thompson) 
regime of scattering the cross-section is large, the energy transferred to the 
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secondary photon is quite small (see Table 3). Consequently, the mean free 
path of these electrons increases with reduction of their energy as 1/e. These 
electrons produce large number of 7-rays with energy below the effective pair- 
production threshold 0, which therefore penetrate very large distances without 
interacting with the ambient background radiation. Consequently, the number 
of these photons remains almost constant. 





Fig. 9. Cascade curves of electrons (upper panel) and photons (bottom panel) for 
showers produced in the Planckian gas of ambient photons by high energy photons. 
The results are obtained for K t h = Sth^o = b and several values of kq = £0^0 
indicated at the curves. 



3 In the case of the black-body radiation of background photons there is no strict 
kinematic threshold; 7-rays with energy e < m e c 2 /3kT may still interact with the 
background photons from the Wien tail of distribution. 
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Fig. 10. Cascade curves of electrons (upper panel) and photons (bottom panel) for 
showers produced in the Planckian gas of ambient photons by high energy photons. 
The results are obtained for kq = eo^o = 10 3 and for values of K t h = e^tio indicated 
at the curves. 

4-1.3 Magnetic field 

Fig. 11 shows the cascade curves of electrons in the magnetic field obtained 
with the adjoint equation technique. The comparison with the Monte Carlo 
simulations [24] shows nice agreement between results obtained by two differ- 
ent methods. 

Following to Ref. [24] we express the depth of penetration of cascade particles 
in units of the radiation length in the magnetic field 

XP = 0.207 • 10- 7 ^Xo /3 , (23) 
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where xo = £oH/H cr . Similar to the cascade in the photon gas, the radiation 
length in the magnetic field depends on the primary energy, but in this case 
the dependence is slower, ~ Eq 3 . 




Fig. 11. Cascade curves of electrons for showers initiated by primary electrons in the 
magnetic field. Different values (indicated at the curves) of the ratio of primary and 
threshold energies are assumed for the fixed xth = £th,H/H CT = 10 3 . For comparison, 
the results obtained in Ref. [24] are also shown (dashed curves). 

Fig. 11 illustrates dependence of cascade curves on the primary energy in a 
deep quantum regime with the threshold value of the parameter Xth = 10 3 . 
It is seen that the location of the cascade curve maximum t m ranges within 
2 3 radiation lengths. In geometrical units this implies that the position of 
the cascade curve maximum increases with energy proportional to e]j 3 . The 
maximum electron number, N m , grows with e as ~ Eq 7 . Thus, the cascade 
development in the magnetic field in the quantum regime is somewhat inter- 
mediate between the cascades in matter and the photon gas. 

In Fig. 12 we present the cascade curves of electrons and photons for different 
values of the threshold energy e th . It is seen that in the regime Xth <^ 1 the 
cascade curves have a behaviour quite similar to that for the cascade curves 
in the photon gas shown in Figure 10. This can be explained by similarities 
of the synchrotron radiation and the Compton scattering in the non-quantum 
regime. In both cases the electrons suffer a large number of collisions but with 
small energy transfer to secondary photons. In addition, in both cases 7-rays 
stop to interact effectively with the ambient medium. 

Table 5 demonstrates the dependence of the N^/N e ratio on xth and t. It 
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Fig. 12. Cascade curves of electrons (upper panel) and photons (bottom panel) for 
showers initiated by primary electrons (solid curves) and photons (dashed curves) 
in the magnetic field. Different values of the threshold energy (eth) are assumed for 
the fixed Eq = 2- 10 and H = 10 11 G. The £o/ £ th ratios are indicated at the curves. 
The corresponding values for xth are presented in Table 5. 

Table 5 

The ratio N 1 /N e for cascade curves presented in Figure 12. 



£o/ £ th 


10 2 


10 4 


10 5 


10 6 


10 8 


Xth 


4.5 • 10 2 


4.5 


4.5 • 10" 1 


4.5 • 10~ 2 


4.5 • 10~ 4 


t =0.3 r.l. 


0.76 


0.92 


1.44 


2.66 


5.95 


t =0.9 r.l. 


0.94 


1.09 


2.73 


4.65 


1.92 


t =3 r.l. 


1.00 


1.17 


11.2 


18.6 


45.2 


t =6 r.l. 


1.03 


1.19 


107 


140 


455 
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behaves quite similar to cascades in matter. Particularly, for large threshold 
energies £ t h the electron and photon numbers are comparable; for small £ t h 
the photon number is considerably larger, and the ratio N y /N e increases with 
the depth. 



4-2 Energy spectra of cascade particles 



Compared to the cascade curves, the energy spectra of electrons and photons 
at different stages of the cascade propagation contain more circumstantial 
information about the shower characteristics. 




Fig. 13. Differential energy spectra of electrons (solid curves) and photons (dashed 
curve) for cascades initiated in the hydrogen gas by a primary photon with energy 
£o = 2 • 10 8 . Different values of depths in units of the cascade maximum t m are 
shown at the curves. 

In Fig. 13 we show the energy spectra of cascade electrons and photons at 
different penetration depths in the hydrogen gas. At energies exceeding the 
critical energy, e ^> £ cr , the spectra of both electrons and photons are de- 
scribed by power-law dN/de ~ e~ a , where the spectral index is function of 
the penetration depth. Near the shower maximum, a ~2 for both electron and 
photon spectra. With depth the spectra become steeper. 

Below the critical energy, the cross-sections of the bremsstrahlung and pair- 
production processes are not sufficiently large to support the cascade devel- 
opment against the dissipative processes like ionization and Compton scatter- 
ing. Thus the multiplication of cascade particles is dramatically reduced. The 
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Fig. 14. Differential energy spectra of cascade electrons (upper panel) and photons 
(bottom panel) in the cascade initiated by a primary photon (kq = 10 3 ) in the 
radiation field with Planckian spectral distribution. The spectra are calculated for 
different penetration depths indicated (in units of radiation lengths) at the curves. 

cooling time of electrons due to ionization losses is proportional to energy, 
therefore below the critical energy this process leads to significant hardening 
of the electron, and correspondingly also the photon spectrum that behaves 
as dN/de oc e^ 1 . 

In Fig. 14 we show the differential energy spectra of cascade particles in the 
radiation field with Planckian type spectral distribution. These spectra are 
quite different from the ones that appear in the cascades developed in matter. 
In the high energy (Klein-Nishina) region, k = eujq ^> 1, and for not very 
large depths (t < 10 r.l.) the shape of the differential spectra of electrons and 
photons is quite insensitive to the depth (contrary to the cascade in matter). At 
these energies the differential spectra are characterized by slopes with indices 
between 1 and 1.5. 
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At k ~ 1 there exists a pronounced transition region where the energy spectra 
undergo dramatic changes. The reason for the appearance of the transition re- 
gion is the change of the character of Compton scattering (from the Thompson 
to the Klein-Nishina regime) and the sharp reduction of the pair production 
cross-section. Obviously, a broader energy distribution of background photons 
would make this transition smoother and less pronounced. 





Fig. 15. Differential energy spectra of cascade electrons (upper panel) and photons 
(bottom panel) in showers initiated in the magnetic field of intensity H = 10 11 G 
by a primary photon of energy eo = 2 • 10 7 . The spectra are calculated at different 
depths indicated (in units of radiation lengths) at the curves. 

At low energies (k ~ 1) the cascade development is not supported by pair- 
production. Here we deal with an ensemble of photons not interacting with the 
environment and an ensemble of electrons continuously cooling down with the 
characteristic (Thompson) time ~ \je. This results in standard electron 
spectra with spectral index a — 2. During propagating into deeper layers of 
the photon gas, these electrons produce large number of low energy photons 
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with energy spectrum dN^/de oc e L5 . 

For any finite depth t, electrons do not have enough time to be cooled down 
to energies e — > 0. Therefore the electron spectrum drops dramatically below 
some energy e e which decreases with depth as ~ I ft. This effect is clearly 
seen Fig. 1 4. The corresponding response in the photon spectrum is also quite 
distinct. Below the break energy around ~ ooo^l, the photon spectrum becomes 
extremely hard, diV 7 /de ~ const. 

The energy spectra of electrons and photons produced during the cascade 
development in the magnetic field are quite similar to the spectra of electro- 
magnetic cascades in the radiation field with Planckian distribution of target 
photons. These spectra are shown in Fig. 15. All features discussed above in 
the context of cascading in the photon gas, are clearly seen also in cascades 
developed in the magnetic field, if we express the penetration depths in units 
of radiation lengths, and the energies of electrons and photons in the form 
of products eu)q and e(H/H cr ). The cascade spectra in radiation and mag- 
netic fields are not, however, identical. For example, because of significant 
differences in the asymptotics of relevant cross-sections, 7-ray spectra in the 
magnetic field in the quantum regime are flatter than the corresponding 7-ray 
spectra in the photon gas in the Klein-Nishina regime (compare Figs. 14 and 
15). 



5 Mixed environment 

In order to reveal peculiarities of the cascade development in different sub- 
stances, in previous sections we limited our discussion to the "clean" environ- 
ments dominated by matter, radiation or magnetic fields. The "pure cascade" 
concept is not only a convenient theoretical approximation. In fact, under cer- 
tain realistic conditions and within limited energy regions, this could be the 
most likely realization of particle interactions with then ambient medium. The 
relativistic electron-photon cascades in the Earth's atmosphere, in the inter- 
galactic medium and in pulsar magnetospheres are 3 characteristic examples 
of "pure" cascade developments in matter, photon gas and magnetic field, 
respectively. 

In some cases, however, "parallel" interactions of electrons and photons with 
2 or 3 substances can proceed simultaneously and with comparable efficien- 
cies. The outcome of the interference of several competing processes could 
be quite different and complex depending on relative densities of the ambi- 
ent plasma, radiation and magnetic fields, as well as on the energy of pri- 
mary particles. For example, interactions of > 10 20 eV protons with 2.7 K 
CMBR leads to production of secondary electrons, positrons and 7-rays of 
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Table 6 

Parameters of the mixed environment used for calculations of the cascade curves 
shown in Fig. 16. H is the strength of the magnetic field, T is the temperature of 
the black-body radiation, wp and wq are energy densities of the magnetic field and 
black-body radiation, respectively. 



Parameter Combination 


1 


2 


3 


4 


5 


H, Gauss 


100 





100 


100 


10 


kT, eV 




3 


0.3 


3 


0.3 


wf, erg/cm 3 


400 





400 


400 


4 


wg, erg/cm 3 





~ 10 4 


~ 1 


~ 10 4 


~ 1 



energy > 10 19 eV, which in their turn trigger electromagnetic cascades in 
the same radiation field. However, due to the synchrotron cooling of elec- 
trons this process would be significantly suppressed if the intergalactic mag- 
netic field exceeds 10~ 10 G. The characteristic energy of synchrotron photons 
~5x 10 8 (S/10 -10 G)(E/10 19 eV) 2 eV is too small for interactions with the 
diffuse photon fields on the Hubble scales. However, if the energy of electrons 
and/or photons injected into the intergalactic medium exceeds > 10 21 eV (this 
could be the case, for example, of secondary products from decays of the so- 
called topological defects), then the synchrotron photons appear in the TeV 
energy range. The TeV 7-rays effectively interact with the diffuse extragalac- 
tic infrared radiation, and trigger new, low energy electron-photon cascades. 
In both cases the synchrotron radiation changes the character of the cascade 
development in the radiation field, but cannot support its "own" cascade in 
the magnetic field. 

Below we briefly discuss a more interesting scenario, when the cascade de- 
velops both in the radiation and magnetic fields. Let's assume that a 7-ray 
photon of energy E y = 10 20 eV is injected into a highly magnetized low-density 
plasma with thermal radiation density comparable to the energy density of the 
magnetic field H 2 /87r. In principle such a situation can occur in the vicinity 
of central engines of AGN. 

Fig. 16 illustrates the impact of the magnetic field and the temperature of 
blackbody radiation on cascade curves. Five different combinations of these pa- 
rameters presented in Table 6 have been analyzed. The minimum particle en- 
ergy for the cascade curves was taken 5 x 10 11 eV, i.e. comparable or larger (for 
the assumed radiation temperatures) to the effective pair production threshold 
in the black-body radiation , ~ m 2 c c 4 /kT ~ 10 n (A;T/l eV)" 1 eV. 

For the chosen parameters, the radiation length in the magnetic field is much 
smaller than the mean free path of primary 7-rays in the photon field. At ab- 
sence of magnetic field (parameter combination 2) primary 7-rays penetrates 
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Fig. 16. The cascade curves of electrons (upper panel) and photons (bottom panel) 
with energy E > E t ^ = 5 • 10 11 eV for showers initiated by a primary photon of 
energy Eq = 10 20 eV in a mixed environment consisting of the magnetic field and 
the blackbody radiation. The calculations are performed for 5 combinations of pa- 
rameters characterizing the target radiation and magnetic fields. These parameters 
are listed in Table 6. 



very deep into the source without interacting with the ambient photon gas. 
Consequently, the cascade starts very late and remains as underdeveloped. 

The presence of a strong magnetic field makes the cascade development much 
more effective, which now is supported by magnetic pair production and syn- 
chrotron radiation. The cascade develops in this regime until the energy of 
7-rays is reduced down to the effective threshold of the magnetic pair pro- 
duction, E$IP ~ 10 17 (#/100 G)" 1 eV. Therefore, in the case of pure magnetic 
field we have quite simple and predictable cascade curves (parameter combi- 
nation 1). 
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The presence of photon gas changes significantly the character of cascade de- 
velopment. At energies above E^ h both electrons and 7-rays interact mainly 
with magnetic field. When the cascade particles enter the energy interval de- 
termined by the pair-production thresholds in the radiation and magnetic 
fields, [Efa \ EtP], the 7-rays interact effectively with the photon gas, while 
the electrons continue to interact mainly with the magnetic field. Although for 
kT = 3 eV, the blackbody radiation density considerably exceeds the energy 
density of 100 G magnetic field (see Table 6), because of the Klein-Nishina 
effect the interactions of electrons with photon gas become significant only 
when the electrons are cooled down to energies < 100 GeV. The mean free 
path of 7-rays in the black-body radiation field has minimum around i? th . 
Therefore, 7-rays start to interact intensively with the ambient photons at 
depths comparable with the mean free path A min ~ A 7 (£' t ^). This results 
in a rapid growth of the electron cascade curve and reduction of the high 
energy photon cascade curve (parameter combinations 3-5). Due to the rapid 
synchrotron cooling, the mean free paths of very energetic electrons in the 
magnetic field are very short, therefore for each given depth the main source 
of these electrons are 7-rays which can penetrate much deeper. On the other 
hand, due to the same synchrotron losses, the electrons cannot support repro- 
duction of very high energy 7-rays. Therefore, the at depths 3> A min both the 
electron and photon curves are suppressed. 

The increase of the temperature increases the energy density of the photon 
gas (oc T 4 ) and makes interactions of cascade particles with radiation more 
intensive. This results in faster absorption of cascade particles (parameter 
combination 4)- The reduction of the magnetic field (parameter combination 5) 
leads to the increase of the magnetic radiation length, as well as to the increase 
of the magnetic pair production threshold. As a result, the cascade develops 
slower. 

Figure 17 illustrates the spectral evolution of cascade 7-rays calculated for 
H = 100 G and two different temperatures of the black-body radiation, 
kT = 0.3 and 3 eV (parameter combinations 3 and 4)- It is seen that at 
small depths (t < 10 10 cm for the upper panel and t < 10 8 cm for the bottom 
panel) the spectra are fully determined by interactions with the magnetic field. 
At these depths a "standard" power-law spectrum with a — 1.5 is formed in 
a broad energy region below E^ h . At larger depths the interactions with the 
photon gas start to deform the shape of the spectrum. These interactions lead 
to absorption of photons above E^\ At the same time, synchrotron radiation 
of the photo-produced pairs appears in the energy region E < Thus, 
only the photons with energy below E^ can survive at large depths. The 
spectrum of these photons is close to power-law with photon index 1.8-1.9, 
i.e. steeper than the the canonical e -15 cascade photon spectrum formed in 
the pure photon gas or pure magnetic field. This is explained by the break 
of symmetry between electrons and 7-rays. While at late stages of the cas- 
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cade development 7-rays continue to produce high energy electron-positron 
pairs, the synchrotron cooling of these electrons does not anymore support the 
cascade, but rather destroys it. 

Finally, in Fig. 18 we show the spectral evolution of cascade 7-rays calculated 
for a less extreme combination of model parameters. Namely, compared to 
Fig. 17, we assume a smaller energy for primary 7-rays (E = 10 17 eV) and 
weaker magnetic field (H = 1 G). At these conditions the primary photons do 
not interact with the magnetic field. Instead, the cascade development starts 
with pair production in the radiation field. Nevertheless, we can see that many 
basic features, in particular the e~ 2 type spectrum at low energies, are quite 
similar to the spectra shown in Fig. 17. 





Fig. 17. Differential energy spectra of 7-rays of the cascade initiated by a primary 
photon of energy Eq = 10 20 eV in the compound environment consisting of a black- 
body radiation of temperature T and a homogeneous magnetic field H=100 G. 
Upper panel - kT = 0.3 eV; bottom panel - kT = 3 eV. The spectra are calculated 
at different depths (in cm) indicated at the curves. 
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Fig. 18. The same as in Fig. 17, but for the energy of primary photons Eq = 10 eV 
and the magnetic field H = 1 G. The temperature of the black-body radiation is 
assumed 0.3 eV. 



6 Summary 



In this study, the technique of adjoint cascade equations has been applied 
to investigate properties of electron-photon cascades in hydrogen gas and in 
ambient radiation and magnetic fields. We also have inspected the main fea- 
tures of cross-sections of relevant processes that initiate and support cascade 
developments in these substances. 



The cascade curves of electrons and photons in the photon gas and magnetic 
field have features quite different from the cascade curves in matter. The 
energy spectra of cascade particles are also considerably different from the 
conventional cascade spectra in matter. The spectra for the magnetic field have 
properties intermediate between those for cascade spectra in matter and in the 
photon field. Although for certain astrophysical scenarios the development of 
cascades in "pure" environments can be considered as an appropriate and fair 
approximation, at some conditions the interference of processes associated 
with interactions of cascade electrons and 7-rays with the ambient photon 
gas and magnetic field (or matter) can significantly change the character of 
cascade development, and consequently the spectra of observed 7-rays. The 
impact is very complex and quite sensitive to the choice of specific parameters. 
Therefore each practical case should be subject to independent studies. 
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